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摘 要 : 有 机 朗 肯 循环 ORC) 技术 可 以 有 效 回收 重型 车 用 柴油 机 的 排 气 热能 。 但 在 道路 运行 过 程 中 排 气 热源 
不 稳定 ， 冷 源 的 维持 也 需要 消耗 能 量 ， 使 得 蒸发 压力 和 冷凝 压力 的 在 线 优 化 与 控制 实施 均 存 在 严峻 挑战 。 为 此 ， 建 立 
了 ORC 余 热 回收 系统 的 效率 预测 模型 ， 以 ORC 系 统 的 热力 学 第 一 定律 模型 为 基础 ， 抓 住 蒸发 器 阶 跃 响应 的 宏观 特 
征 ， 提 出 了 有 效 换 热 量 的 动态 修正 模型 ， 以 及 工 况 模式 切换 的 预 判 因子 ， 并 采用 简化 机 理 和 数据 建 模 相 结合 的 方法 ， 
建立 了 冷凝 器 散热 功 耗 的 经 验 模型 ， 提 出 了 基于 效率 预测 模型 的 ORC 余 热 回 收 系统 优化 控制 架构 及 策略 。 道 路 仿真 
结果 表明 了 效率 模型 和 优化 算法 的 有 效 性 : ORC 系 统 的 有 效 做 功 时 间 达 到 94%， 蒸 汽 过 热度 控制 在 5-15K 之 间 ， 冷 凝 
器 散热 功 耗 与 发 电功率 的 比值 维持 在 20% 以 下 。 
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Efficiency Prediction Model of a Diesel-Engine ORC Exhaust Heat Recovery 


System 


YANG Can, XIE Hui, ZHOU Kui, CHENG Zhang 
(State Key Laboratory of Engines, Tianjin University, Tianjin 300072, China) 


Abstract: The ORC (organic Rankine cycle) System is regarded as one of the most potential ways for 
exhaust heat recovery of heavy-duty vehicle diesel engines. But the exhaust energy fluctuates widely during 
driving cycles, also the ORC cooling will introduce extra energy consumption. Consequently, it is a big 
challenge to optimize and control the evaporation pressure and condensation pressure in the driving cycle. 
To this end, a control-oriented ORC efficiency prediction model was completed at first: basing on a 
foundation model of the first law of thermodynamics; focusing on the dominant feature of the ORC system 
step response, an effective heat absorption dynamic model and operating modes corrected factor were built 
up; and an empirical model was established to describe the energy consumption for ORC cooling. 
Afterwards, a sequence optimization algorithm was proposed to seek optimum values of the evaporation 
pressure and condensation pressure dynamically. The validation results based on a simulation plat were 
encouraging: the effective working time percentage achieves 94%, the superheat degree is controlled as SK- 
1SK, and the cooling power consumption 1s kept less than 20% of the ORC power generation. 
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0 引言 


耗 法 规 的 热点 技术 m]。BMW、AVL、 橡 树 岭 等 研 
完结 果 表 明 ， 该 技术 在 重型 车 上 可 获得 5-10% 的 
燃油 经 济 性 改善 效果 中。 


采用 ORC 技术 回收 和 利用 重型 柴油 机 的 余 寺 


注 


能 是 改善 发 动机 燃油 经 济 性 ， 满 足 未 来 更 严 苛 ; 但 车 辆 在 道路 运行 过 程 中 ， 其 发 动机 基本 处 
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es 大 负荷 、 小 负荷 等 动态 变 工 况 运行 状 

， 内 燃 机 的 排 气 流量 和 温度 也 必然 呈现 出 强烈 
re 前 期 研究 67 发 
现 : ORC 余热 回收 系统 在 循环 工 况 中 呈现 出 强烈 
的 动态 行为 特征 ， 包 括 蒸 汽 参数 的 波动 以 及 工 况 
模式 的 切换 ， 从 而 恶化 了 其 实际 的 节 油 效果 。 法 
国 石油 研究 院 IFP 的 Tonalsl 等 人 研究 指出 ORC 系 
统 在 动态 变 工 况 条 件 下 的 优化 控制 是 其 应 用 的 主 
要 瓶颈 。 

此 外 ， 热 机 不 可 能 从 单一 热源 取 热 ， 并 使 之 
完全 转化 为 有 用 功 而 不 产生 其 它 影响 。ORC 系统 
在 回收 发 动机 排 气 能 量 发 电 的 同时 ， 也 通过 冷凝 
器 将 低温 热量 释放 到 发 动机 冷却 系统 当中 ， 增 加 
了 冷却 系统 的 散热 负荷 和 散热 功 耗 。Horst TA 等 
中 的 研究 表明 ，ORC 系统 散热 会 增加 冷却 水 泵 的 
功 耗 ， 从 而 导致 其 节 油 潜力 下 降 了 1.3%-1.8%。 
ORC 系统 回收 的 余热 能 越 多 ， 其 散热 需求 也 越 
大 ， 冷 却 系 统 的 散热 功 耗 也 会 越 大 。ORC 系统 发 
电 和 发 动机 冷却 系统 耗 电 两 者 之 间 存 在 着 一 定 的 
三 亏 夺 合 关系 [10]。 

可 以 发 现 ，ORC 系统 在 回收 车 用 内 燃 机 排 气 
能 量 时 ， 会 呈现 出 强烈 的 动态 特征 和 耦合 特征 ， 
使 得 ORC 系统 在 动态 变 工 况 条 件 下 的 优化 决策 和 
控制 实施 均 存 在 严峻 挑战 。 为 此 ， 本 文 以 ORC 系 
统 的 动态 特征 和 耦合 特征 为 核心 ， 则 在 建立 面向 
ORC 系统 优化 控制 的 效率 预测 模型 ， 支 持 对 蒸发 
压力 和 冷凝 压力 进行 在 线 决 策 和 控制 。 


1 问题 描述 
1.1 研究 对 象 
图 1 为 本 文 研究 对 象 的 系统 构 型 : 采用 ORC 


技术 回收 一 台 重 型 车 用 柴油 机 增 压 涡轮 后 的 排 气 
能 量 ;， 并 为 ORC 系统 匹配 了 一 套 可 控 的 冷却 系 


迎风 冷却 


SE, A | 


膨胀 阀 ” 透 平 膨胀 机 


工 质 泵 储 液 镀 冷却 水 泵 


pa 


图 1 柴油 机 ORC 余热 回收 系统 的 构 型 示意 
Fig.1 Schematic diagram of the diesel-engine ORC exhaust 
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heat recovery system 
1.2 问题 描述 
ORC 系统 优化 和 控制 研究 的 核心 任务 是 : 在 
道路 运行 过 程 中 ， 根 据 排 气流 量 和 排 气 温度 的 瞬 
变 特征 ， 以 系统 效率 最 优 为 目标 ， 实 时 决策 ORC 
系统 的 运行 状态 ， 并 通过 控制 手段 加 以 实现 。 
ORC 系统 的 系统 效率 为 
Wc W. 
Ez (1) 
= / (Xic | Rs,) 
式 中 ，www 为 系统 效率 ，E。 为 排 气 能 量 ，WW。 
为 ORC 系统 的 循环 净 功 ，W. 为 冷却 附件 的 功率 
消耗 ， Rs, 表示 排 气 边界 条 件 ， Xi 表示 ORC 系 
统 的 运行 状态 ， 包 括 营 发 压力 、 过 热度 、 冷 凝 压 
力 和 过 冷 度 
要 在 整个 动态 运行 历程 中 ， 时 刻 保持 ORC 
es a *， 获 得 最 优 的 系统 效率 ， 则 
要 在 每 一 个 控制 步 长 ， 根 据 排 气 边界 条 件 ， 包 
括 排 气温 度 和 排 气流 量 ， 优 化 和 控制 ORC 系统 的 
蒸发 压力 、 过 热度 、 冷 凝 压力 和 过 冷 度 。 其 中 ， 
根据 近期 的 研究 @。 2 区 可 以 发 现 ， 过 热度 和 过 冷 度 
的 最 优 取 值 比较 简单 ， 不 会 随 着 运行 工 况 的 变化 
而 变化 。 然 而 ， 营 发 压力 和 冷凝 压力 的 最 优化 规 
律 尚 不 清晰 ， 分 而 论 之 ， 主 要 有 如 下 三 点 原因 : 


Dsystem = 


1.2.1 动态 es 

稳 态 研究 [9 表明 ， 最 优 的 蒸发 压力 与 排 气温 
度 正 相关 。 排 气温 度 和 蒸发 压力 在 数值 上 是 一 一 
对 应 的 ， 但 在 动态 运行 过 程 中 ， 两 者 在 时 间 上 并 
不 一 致 ， 存 在 一 定 的 时 间 差 ， 因 为 发 动机 和 ORC 
两 个 系统 转化 能 量 的 时 间 尺 度 存在 较 大 差异 ， 如 
表 1 所 示 。 


表 1 发 动机 与 ORC 系统 能 量 转化 特征 的 对 比 


Tablel Comparison of energy conversion features between 


engine and ORC 
发 动机 ORC 系统 
装置 工作 方式 往复 连续 
能 量 转化 方式 缸 内 燃烧 间 壁 传 热 
动态 时 间 尺 度 快 : 工作 循环 慢 : 分 钟 


发 动机 往复 式 的 工作 方式 ， 以 及 爆炸 燃烧 的 
能 量 释 放 形 式 ， 使 得 燃油 的 化 学 能 在 一 个 工作 循 
环 内 完成 到 机 械 能 的 转化 ， 其 能 量 转化 的 动态 时 
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间 尺 度 较 快 ， 为 一 个 工作 循环 。 然 而 ，ORC 的 工 
作 方 式 是 连续 的 ， 发 动机 瞬 态 的 排 气 能 量 经 过 间 
壁 传 热 被 ORC 工 质 吸收 ， 加 热 燕 发 ， 推 动 膨胀 机 
做 功 产 生 电 能 ， 系 统 的 传 热 和 流动 惯性 使 得 ORC 
能 量 转化 的 动态 时 间 尺 度 较 大 ， 需 要 几 十 秒 甚 至 
数 分 钟 ， 如 表 1 所 示 。 这 种 动态 特征 使 得 蒸发 压 
力 和 排 气温 度 一 一 对 应 的 稳 态 控制 律 在 动态 运行 
过 程 中 无 法 实施 ， 需 要 对 蒸发 压力 进行 在 线 优 
化 。 


1.2.1 工 况 特性 

Tonals 等 人 研究 指出 ，ORC 系统 从 冷 起 动 开 
始 ， 逐 渐 被 加 热 、 升 压 ， 直 至 完全 蒸发 、 过 热 ， 
到 推动 膨胀 机 对 外 输出 做 功 ， 会 经 历 一 系列 的 动 
态 过 程 。Xielg] 等 人 将 其 定义 为 四 种 基本 的 工 况 模 
式 ， 包 括 起 动 模式 、 冲 转 模 式 、 做 功 模 式 和 保护 
模式 。 并 指出 工 况 模式 的 波动 会 降低 ORC 系统 的 
有 效 做 功 时 间 ， 进 而 影响 节 油 效果 的 发 挥 。 因 
此 ， 在 道路 运行 过 程 中 ，ORC 系统 优化 与 控制 更 
大 的 任务 是 调节 花 发 压力 ， 避 免 ORC 系统 落 入 到 
保护 模式 ， 绥 解 工 况 模式 的 波动 ， 提 高 ORC 系统 
有 效 做 功 时 间 。 因 此 ， 调 节 蒸 发 压力 不 单单 是 优 
化 ORC 系统 在 每 一 个 控制 步 长 上 的 系统 效率 ， 更 
重要 的 是 ， 避 免 ORC 系统 工 况 模式 的 波动 。 


节 汽 过 热度 线 


透 平 转速 线 


包 三 

0 
如 透 平 高 效 区 | 红 
吕 转速 低 阔 值 | 博 
外 1 控 
扩 证 

网 


运行 时 间 [S] 


图 2 ORC 系统 动态 的 工作 过 程 及 模式 定义 
Fig.2 ”Dynamic working process and mode definition of the 
ORC system 


1.2.1 看 合 特性 

热机 不 可 能 从 单一 热源 取 热 ， 并 使 之 完全 转 
化 为 有 用 功 而 不 产生 其 它 影 响 。ORC 系统 在 回收 
发 动机 排 气 能 量 发 电 的 同时 ， 也 通过 冷凝 器 将 低 
温 热量 释放 到 发 动机 冷却 系统 当中 ， 增 加 了 冷却 
系统 的 散热 负荷 和 散热 功 耗 。Yangllo 等 人 研究 表 
明 ，ORC 系统 回收 的 余热 能 越 多 ， 其 散热 需求 也 
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越 大 ， 冷 却 系统 的 散热 功 耗 也 会 越 大 。ORC 系统 
发 电 和 发 动机 冷却 系统 耗 电 两 者 之 间 存 在 强烈 的 
耦合 关系 ， 如 图 3 所 示 。 而 冷凝 压力 的 取 值 对 这 
种 耦合 关系 有 着 强烈 的 调控 作用 ， 应 随 着 发 动机 
转速 和 负荷 的 增加 而 增加 。 换 句 话 说， 冷凝 压力 
的 最 优 取 值 也 受到 排 气 条 件 的 强烈 影响 。 


9 一 ORC 发 电功率 外 包 络 线 
12 -ec- 冷凝 压力 3bar 
十 0 一 冷凝 压力 4bar 
10 一 一 全 冷凝 压力 5bar 
十 一 冷凝 压力 6bar 
8 十 < 一 冷凝 压力 7bar 7 
二 <- 冷凝 压力 8bar A 


冷却 系统 的 功率 消耗 [kW] 


oo 一 
a 
EE 


图 3 冷凝 压力 对 耦合 关系 的 调控 作用 


g.3 Effects ofthe condensing pressure on the profit-loss 


Fi 


上 


relationship 


综 上 所 述 ， 在 道路 运行 过 程 中 ， 因 为 ORC 系 
统 的 动态 特性 、 工 况 特征 和 耦合 特性 ， 使 得 蒸发 
压力 和 冷凝 压力 的 优化 和 控制 存在 严峻 挑战 ， 而 
建立 面向 ORC 系统 优化 控制 的 效率 预测 模型 是 突 
破 这 一 挑战 的 关键 。 


2 效率 预测 模型 


2.1 控制 架构 

综合 考虑 ORC 系统 的 动态 特性 、 工 况 特征 和 
耦合 特性 ， 提 出 了 ORC 系统 的 优化 控制 架构 ， 如 
图 4 所 示 。 针 对 给 定 排 气 边界 条 件 R, ， 管 理 器 和 
效率 观测 器 协同 工作 ， 优 化 得 到 最 优 的 蒸发 压力 
Pe 和 冷凝 压力 必 .,， 和 蒸汽 温度 Te, 和 透 平 转 
速 Wec7y 一 起 作为 控制 器 的 控制 目标 及 w ， 在 控 
制 层 面 直 接 采 用 成 熟 的 PID 控制 算法 ， 控 制 ORC 
系统 透 平 转速 、 膨 胀 比 ， 工 质 泵 流量 ， 以 及 冷却 
系统 附件 功率 等 。 该 架构 的 关键 之 处 在 于 建立 一 
个 宏观 的 、 具 有 预测 功能 的 效率 模型 。 


2.2 效率 模型 

根据 图 4， 效 率 模型 的 核心 功能 是 ， 在 管理 器 
的 任意 决策 时 刻 x， 根 据 该 时 刻 的 排 气 参数 
Rs,(K) 以 及 ORC 系统 的 状态 参数 Xi(x) ， 需 要 预 
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测 ， 如 果 管 理 器 给 出 决策 变量 为 Xic o, (x)， 则 在 
x+1 时 刻 ORC 系统 的 工 况 模式 、 输 出 功率 以 及 冷 


却 功 耗 ， 将 系统 效率 ,w(x+1) 反馈 给 管理 器 。 
系统 效率 15 的 定义 式 为 公式 (2)。 


TRC,T 


管理 器 
工 况 调度 机 


工 况 

效率 优化 器 

效率 模型 
模式 预测 || 功率 预测 | 功 耗 预测 
图 4 ORC 余热 回收 系统 的 优化 控制 架构 
Fig.4 Control framework of the ORC exhaust heat recovery System 
Wh 吸 热 量 的 实时 预测 模型 ， 考 虑 管理 器 的 决策 变量 

sytem = 以 及 控制 器 的 执行 能 力 ， 基 于 ORC 系统 的 热力 学 


Ex (2) 
一 f (Rg,, XRc, Opt» Yic) 


热度 ， 并 根据 冷却 功 耗 模 


模型 ， 预 测 ORC 系统 的 输出 功率 、 


散热 量 以 及 过 


型 和 工 况 模 型 ， 预 测 冷 
却 功 耗 和 工 况 模式 ， 在 此 基础 上 计算 ORC 系统 的 


其 中 ， Bi 表示 发 动机 的 排 洁 能量，Wro(r) 表 不 效率 。 在 效率 模型 的 建立 过 程 中 ， 首 先 作 如 下 两 
ORC 系统 未 来 一 段 时 间 内 的 发 电量 ，.(7) 表示 点 假设 : 
冷却 系统 所 消耗 的 功率 ，&, 表示 ORC 系统 的 工 1 管理 器 给 出 的 决策 变量 ，， 被 控制 器 
人 完全 执行 ， 不 考虑 执行 器 的 执行 过 程 和 
根据 ORC 系统 的 动态 特性 ， 发 动机 工作 循环 执行 能 力 。 
是 其 完成 能 量 转 化 的 基本 时 间 单位 ， 如 果 以 发 动 2) 不 考虑 膨胀 机 、 发 电机 、 工 质 泵 、 冷 却 
机 的 工作 循环 为 时 间 基 准 的 话 ，ORC 系统 的 动态 系统 和 阀门 的 硬件 执行 过 程 ， 只 关注 其 
es be 输出 量 ， 包 括 膨胀 比 x,。, 、 发 电 负载 
个 发 动机 工作 循环 。 因 此 ORC 系统 的 发 电量 un tne 
a a a Tre 、 工 质 流量 mc 、 冷却 附件 功率 
相应 的 ，ORC 系统 发 电 所 需要 的 冷却 功 耗 用 所 和 阀门 开 度 Valves 。 
Wet) 家 二 管理 器 
除 此 之 外 ，ORC 系统 的 效率 模型 还 需要 具备 pev PT 
一 项 重要 的 功能 ， 就 是 预测 ORC 系统 的 工 况 模 — 人 | 
式 ， 其 对 于 系统 效率 1, 而 言 是 一 个 0-1 关系 。 台 ||[ 租 Ee 和 
Ex QF: 


、 ORC 系统 的 工 况 模式 ， 是 避免 其 模式 
波动 ， 提 高 做 功 模 式 所 占 时 间 比 例 的 关键 。 


总 之 ，ORC 系统 的 效率 模型 描述 的 对 象 包括 
ORC 系统 及 其 控制 器 ， 需 要 能 够 动态 地 预测 ORC 


系统 的 工 况 模式 所 、 发 电功率 (7) 以 及 冷却 功 效率 计算 
耗 W.(r)， 其 建 模 思路 如 图 5 所 示 。 ^ 


利用 蒜 发 器 的 传 热 动态 效应 ， 建 立 ORC 系统 


Wrc(D -WT) . 


E, 


© 


5 ORC 系统 效率 预测 模型 的 建 模 思路 
Fig.5 Structure diagram of the ORC system efficiency 


modeling 


2.2.1 有 效 吸 热量 

ORC 系统 的 有 效 吸 热量 受到 蒸发 器 对 流 换 热 
过 程 的 控制 ， 考 虑 到 所 建 效 率 模 型 需要 在 控制 器 
中 实时 计算 ， 因 此 ， 其 模型 不 能 够 过 于 复杂 ， 需 
要 在 描述 精度 和 计算 实时 性 之 间 进 行 折 中 处 理 ， 
然后 将 模型 中 不 准确 的 部 分 用 先进 的 控制 算法 进 
行 补偿 。 因 此 ， 蒸 发 器 对 流 换 热 过 程 的 描述 采用 
基于 集 总 参数 法 的 平均 值 模型 。 
在 排 气 侧 ， 排 气 外 掠 又 排 管束 换 热 过 程 的 控 
制 方程 包括 ， 对 流 换 热 方程 、 能 量 守 恒 方程 和 对 
数 温差 方程 ， 如 公式 (2)、(3) 和 (4)， 


OF 部 hs As, ATr. ( 2 ) 
OF i CEx Ms, (Ts in Ts ee ( 3 ) 
AT an 一 AT mow 
AT 六 = (4) 
. 1 AT Erin 
nn 一 
AT Erou 


其 中 ， 乓 为 排 气 侧 的 平均 换 热 系统 ， 其 服从 茹 卡 
马 斯 卡 斯 (Zhukauskas) 换 热 准则 ， 在 此 基础 上 
可 以 增加 排 气 换 热 系数 的 标定 因子 4. ;4s 表示 
排 气 侧 的 有 效 换 热 面积 ， 罗 mx 表示 排 气 的 质量 流 
量 ， 了 也 


及 和 了 分 别 为 排 气 在 蒸发 器 入 口 处 、 出 
口 处 与 管 壁 平 均 温 度 Tez 的 差 值 ，cs 为 某 油 机 排 


气 的 比热容 ， 工 程 上 常 取 值 为 1.088 kJ/kgK)。 联 
立 求解 可 以 得 到 ， 


_hErAE: 


OF 3 Cp py (T pein 一 Tgp)ll —e ca ) (5$ ) 
公式 (5) 可 以 解读 为 两 个 部 分 ， 最 大 换 热 量 
Op Max 和 排 气 侧 换 热效率 NE E 9 分 别 可 表示 为 
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个 无 量 纲 书 。st 越 大 ， 发 生 于 排 气 与 蒸发 器 管 辟 
之 间 的 对 流 换 热 过 程 就 愈 强烈 。“ 是 排 气 的 换 热 
面积 与 通 流 面积 之 比 。“ 越 大 ， 排 气 与 蒸发 器 管 
壁 之 间 的 换 热 过 程 进行 得 越 充分 ， 换 热效率 越 
高 。 需 要 说 明 的 是 ， 燕 发 器 管 壁 的 平均 温度 7 不 
能 直接 测量 。 但 是 ， 因 为 管 壁 蓄 热 作用 的 存在 ， 
其 温度 不 可 能 快速 变换 ， 因 此 可 以 选 代 求解 ， 在 
x+1 时 刻 ， 可 以 采用 第 x 时 刻 的 管 壁 温度 T(r) 来 
计算 排 气 侧 的 换 热 量 0,(x) ， 然 后 ， 根 据 管 壁 的 
能 量 守恒 重新 计算 x+1 时 刻 对 应 的 管 壁 温度 
Tip(K+1)。 

在 计算 得 到 排 气 换 热 量 Q， 之后， 可 以 根据 燕 
发 器 管 壁 的 集 总 传 热 模型 ， 计 算 ORC 工 质 的 有 效 
吸 热量 On ,; ， 其 控制 方程 为 


By dTgp _ 
Epigp 村 
dT 


QF.— Qkge., eff (9) 


Qnc wy = hrc asArc a(Tap—Trcr) (10) 

如 前 面 所 述 ， 在 控制 器 功能 执行 完美 的 假设 下 ， 
液 相 、 两 相 和 汽 相 的 边界 相对 稳定 ， 故 可 认为 工 
质 在 蒸发 器 中 换 热 的 平均 温度 nc 受到 蒸发 压力 
Pi ;唯一 控制 ， 蒸 发 压力 高 ， 则 平均 温度 高 ， 
Tac.z 是 关于 决策 量 PB ;的 函数 。 在 这 样 的 条 件 
下 ， 公 式 (9) 和 公式 (10) 无 法 得 到 解析 解 ， 微 
分 方程 组 联 立 求解 的 计算 负担 是 内 燃 机 控制 器 无 
法 承担 的 ， 必 须 得 换 一 种 角度 来 解决 这 个 问题 。 
假设 第 * 时刻 的 蔡 发 压力 为 Pe 3(x) ， 在 x+l 
时 刻 管理 器 给 出 蒸发 压力 决策 值 为 Pu :(c+D ， 并 
且 Pic s(x+])zPic s(x)， 那 么 工 质 的 平均 温度 必然 
会 发 生 相 应 的 变化 ， 使 得 Txc a(x+D) Txc. eR 。 
然而 ， 因 为 蒸发 器 管 壁 的 蕾 热 作 用 ， 其 温度 To 不 
会 立即 发 生 改变 。 根 据 〈10) 可知， 其 有 效 吸 热 
量 Ou wwr(k+D 会 因为 工 质 平均 温度 Tacs(k+D 的 
突变 而 产生 一 个 脉冲 响应 ， 脉 冲 能 量 的 大 小 等 于 


OR = Cp Mg, (T gm— TP) (6) 
Nag =1-e oi (7) 


在 公式 (7) 中， 指数 项 表征 了 排 气 的 实际 换 热 热 


蒸发 器 管 壁 因 为 温度 变化 1 
量 。 如 果 Pcs3(K+l)< Pc s( 


j 释 放 或 者 吸收 的 能 
x) ， 则 蒸发 器 管 壁 释放 


能 量 ， 反 之 则 吸收 能 量 ， 


其 值 表 征 为 


流 与 可 传递 的 最 大 热流 密度 的 比值 ， 其 可 进一步 
表示 为 ， 


nee=l-e (8) 


其 中 ，& 为 斯 坦 顿 数 ， 是 表征 强制 对 流 强 弱 的 一 


TD 


OO 网 = CgpmgpATEep (11) 
其 中 ，ATss 可 以 表示 为 蒸发 压力 Pi ; 变化 而 导致 
工 质 饱和 温度 变化 量 AT ,的 线性 关系 ， 系 数 为 
Ap 。 通 过 这 种 处 理 ， 可 将 公式 〈9) 和 (10) 改 


_ Qa+O(Qgp) 

Orc, 六 l+T cgr 
其 中 ，7w 表示 蒸发 器 动态 传 热 的 惯性 时 间 常 
数 ， 表 征 了 蒸发 器 的 蕾 热能 力 0 
比值 。 考 虑 到 工 质 侧 换 热 系数 预测 的 复杂 性 ， 其 
惯性 时 间 7ex 可 通过 激励 试验 得 到 ， 并 且 考 虑 到 
流量 和 温度 对 雷诺 数 和 物性 参数 的 影响 ， 故 对 激 
励 试验 得 到 的 惯性 时 间 进 行 流 量 和 温度 修正 ， 修 
正 因子 为 Xrc 和 Xrcr。 在 此 基础 上 ， 根 据 公式 
(5) 和 公式 〈12) 可 以 实时 预测 ORC 系统 燕 发 
器 的 有 效 吸 热量 Ojc ur(z) 。 


(12) 


2.2.2 工 况 模式 因子 
工 况 模式 上, 是 一 个 0-1 变量 ， 在 有 效 吸 热量 
Qxcw 给 定 的 情况 下 ， 受 到 决策 变量 Pac: ， 以 及 
di ee a 多 响 ， 但 是 管理 
需要 知道 当 决 策 变量 Pac.; 被 执行 ， 且 质量 流量 
He ne 
是 否 有 可 能 沙 入 到 保护 模式 ， 关 注 的 是 一 种 可 能 


对 于 给 定 的 有 效 吸 热量 Ouc jy ， 以 及 决策 变 


量 Prc,3， 计 算 当 决策 变量 Pac. 和 控制 器 调节 质量 


流量 mx ,来 不 及 影响 蒸汽 状态 时 ，ORC 系统 瞬时 
的 工 况 模式 上 。 也 就 是 说 ， 在 第 < 时 刻 ，ORC 系 
统 工 质 在 蒸发 器 入 口 的 温度 Tac (x) 和 质量 流量 
mnc r(G) 与 第 x-1 时 刻 的 值 相同 ， 那 么 要 使 蒸发 
器 出 口 薰 汽 的 温度 达到 营 发 压力 Prc,; 对 应 的 安全 
温度 Tkc ;ws， 则 需要 的 最 小 能 量 Oc oy, (x) 为 


Onc ef min(K) = Mage, V (x D(hrec. 3min 


则 工 况 模式 上 为 


hac s(x —1)) 


GlK) 县 if (Qke. eff (x) 本 Onc e 六 mi (x)) (14) 
其 含义 为 ， 如 果 有 效 吸 热量 Oy(x) 大 于 
Qrc wmn(x) ， 则 ORC 系统 工作 在 做 功 模式 ， < 


取 值 为 1， 反之 ORC 系统 会 进入 到 保护 模式 ，< 
取 值 为 0。 


2.2.3 冷却 消耗 功率 
ORC 系统 发 电功率 历 rc(7) 表示 当 决 策 变量 
Pac,3 和 PRc,4 在 一 个 决策 周期 内 的 被 控制 器 完全 执 
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行 之 后 的 预测 值 。ORC 系统 中 工 质 泵 入 口 为 状态 
1 点 ， 如 图 6 所 示 ， 其 饱和 压力 为 Pacy ， 其 比 烩 
为 瑟 kci; 根据 工 质 泵 的 效率 wc ps 可 以 计算 得 到 
2 的 状态 ， 用 一 对 变量 ( prc 3，Trcs ) 表示 ， 其 


比 灼 为 太 ic.;，; 3 点 的 状态 由 决策 变量 pic 控制， 
一 对 变量 ( Pre3» Tac,3); 其 比 烩 为 豆 xc 3; 最 


后 ， 根 据 膨 胀 机 的 等 炉 效率 wi ys 可 以 计算 得 到 4 
点 比 始 为 且 key， 在 此 基础 上 ， 调 用 REFPROP 软 


件 ， 可 以 计算 得 到 ORC 系统 在 该 条 件 下 的 循环 效 
率 为 


(Hie,s 于 


Hic.a)—(Haes —Hia) 
H -HH 


Nrc , Cycle 


MM RC2 


则 ORC 系统 发 电功率 的 预测 值 为 
Wi (z) = Qkc. ef (Dre Cycle (16) 
同时 ， 冷 凝 器 的 散热 量 为 


Qe.con(T) = Onc of (7T)(1 -rc ce) (17) 


6 ORC 系统 循环 示意 图 
Fig.6 Thermodynamic cycle diagram of the ORC system 


el 念 却 消耗 功率 主要 包括 水 和 泵 功 耗 、 
II 
对 于 水 泵 ， 其 功 耗 可 表示 为 
W,, = AV,AP (18) 


式 中 ，%1 是 水 有 泵 的 效率 的 倒数， ee 
水 泵 的 体积 流 率 ，AP,, 表示 水 泵 的 实际 扬程 。 当 
然 ， 对 于 冷却 水 的 流动 循环 而 言 ， 水 泵 的 扬程 其 
实 是 为 了 克服 管 路 阻力 ， 包 括 沿 程 损失 和 局 部 损 
失 。 因 此 ， 水 泵 的 扬程 AP 可 近似 表示 为 
AFp = stp (19) 
将 忽略 一 次 项 所 产生 的 误差 放 在 系数 即 中 进行 标 
定 ， 故 整合 公式 〈18) 和 公式 〈19) 可 以 得 到 
W., = Fo (20) 


(15) 


一 般 地 ， 冷 却 水 流量 和 散热 量 直接 相关 ， 可 以 将 
公式 〈20) 变形 为 
W,, = A Oc, (21) 

式 中 ，0, 为 冷凝 器 中 的 散热 量 。 在 公式 (21) 
的 基础 上 ， 可 以 将 、 罗 和 妙 合 并 ， 得 到 
Wo, = hcp Oe (22) 

同 理 ， 可 以 推导 得 到 风 肩 功 耗 ， 以 及 迎风 阻 
力 所 产 生 功 耗 的 表征 模型 ， 如 公式 (23) 和 
(24)， 


Wo = LO Our (23) 


Wo, = Wes On s (24) 
其 中 ，0O% ;和 Qn 。 分别 表 示 风 局 冷却 和 迎风 冷 
却 从 散热 水 箱 中 带 走 的 热量 。 迎 风 冷 却 是 被 动作 
用 的 ， 受 道路 工 况 影响 ， 而 风扇 辅助 冷却 是 可 以 
主动 调节 的 ， 但 两 者 的 效率 相差 很 大 。 如 果 是 风 
扇 辅 助 冷 却 ， 那 么 其 功 耗 为 


Wr = 至 BV.sAP, (25) 
如 果 是 迎风 冷却 ， 其 阻力 所 产生 的 功 耗 近似 为 
W,, = 记 , AP, (26) 


因为 迎风 冷却 是 直接 作用 ， 没 有 能 量 转 化 环节 ， 
效率 近似 为 100%; 而 风扇 是 一 个 效率 较 低 部 件 ， 
其 等 粹 效率 通常 在 30% 以 下 。 引 入 一 个 迎风 冷却 
的 比例 因子 w, ， 表 征 了 ORC 系统 与 发 动机 冷却 
系统 的 集成 程度 。 当 迎风 冷却 所 占 比例 增高 时 ， 
ORC 放 热 过 程 的 能 效 比 增 大 。 将 公式 (23) 和 
(24) 合并 ， 得 到 
W. = f(@)Q% (27) 
根据 能 量 守恒 ， 冷 凝 器 的 散热 量 和 散热 水 箱 
的 散热 量 相等 ， 即 
Onc. eon = Qo0n = Q paa (28) 
根据 公式 〈22) 和 “(27)， 可 以 得 到 放 热 过 程 
总 功 耗 的 特征 模型 
W. = (hc + f (Or)Qre mw (29) 
定义 公式 〈29) 中 子 项 “+f(@,) ”的 倒数 
为 ORC 冷却 系统 的 能 效 系 数 COP,,, ， 受 冷凝 压 
力 控制 。 根 据 Yangt19 等 人 的 研究 结果 ， 当 全 部 采 
用 风扇 冷却 时 ， 
COP,,, =-3675x phc s+68420x pacs 
-285200x P_ ,+429300 -0 
2.2.4 效率 模型 集成 
根据 效率 模型 的 建立 过 程 ， 将 公式 (12) 采 
用 后 向 离散 处 理 可 以 得 到 
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1 J 


Onc,op (K)—(1— Te Fe op (K— D 


(31) 


= (Qi.(K)+ 6(Qs, (x) 


其 中 ，Az 为 决策 步 长 ， 第 x 时 刻 排 气 的 放 热 量 
Qa.(K) 为 

OQ (K) = cm 
式 中 ， 蒸 发 器 管 壁 的 平均 温度 采用 第 <-1 时 刻 的 


温度 值 ， 而 第 x 时 刻 的 温度 值 可 由 公式 〈33) 进 
行 计算 ， 
CgpMgp (T gp(K) 到 T gp(K 3 D) 
(33) 
= OF, (K)— Once. ef (x) 
第 <x 时刻， 因为 决策 变量 Pc ,(x) 变化 导致 蒸发 器 


管 壁 热 容 吸 收 或 者 释放 的 脉冲 能 量 56(Ou(x)) 为 
O(Qgp(K)) = CgpmgphepAT RC,S (x) (34) 

式 中 ，ATic s(x) 表示 第 * 时 刻 决策 变量 P(x) 对 
应 的 工 质 的 饱和 温度 ， 与 第 x-1 时 刻 决策 变量 
Pei(k-D 对 应 的 工 质 的 饱和 温度 之 差 。 

并 根据 〈14)、(16) 和 〈29) 进行 分 别 计算 
第 < 时 刻 的 工 况 模式 < (x) 、 发 电功率 球 xc(c) 和 冷 
却 功 耗 矿 c(c) ， 其 离散 形式 为 
GK)=if (Orc og (K)— Qc op min(K)) (35) 


Wi (K)= Qc., ef (Kc Cycle (x) (36) 
W(x) = Ee of (zx)(1— Nec, wae (#)) | (37) 
COF,,, (x) 
由 此 ， 可 计算 ORC 系统 在 第 x 时 刻 的 系统 效率 为 
四 Wic (Kx)— W(x) 
WN sysiem — ”加 Ce E,(K) (38) 


总 共有 四 个 系数 需要 标定 ， 分 别 为 排 气 侧 换 热 系 
数 缩放 因子 2 总， 惯性 时 间 7om 的 流量 修正 因子 
Xrcwn 和 温度 修正 因子 Xrc,r， 以 及 蒸发 右 管 壁 热 
容 吸 收 或 者 释放 的 脉冲 能 量 温差 修 正 因子 4 。 


3 效率 模型 验证 


基于 文献 I 所 示 仿 真 平台 ， 对 所 建 效率 模型 
进行 验证 ， 并 且 应 用 到 图 4 所 示 的 控制 架构 中 ， 
对 ORC 系统 进行 在 线 优化 和 管理 ， 从 应 用 效果 上 
进一步 验证 所 建立 效率 模型 的 有 效 性 。 


3.1 模型 精度 验证 
在 低 转 速 小 负荷 ， 到 高 转速 大 负荷 四 个 工 况 
点 对 效率 模型 进行 标定 和 校准 ， 参 数 设置 如 表 2 


Mg (K)(T goin(K) es T gp(K 地 1D))7 ge,. E ( 32) 


所 示 。 效 率 模 型 的 校准 结果 如 图 7 所 示 。 


冷却 方式 


ORC 统 效率 [%] 
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表 2 参数 设置 
Table2 Parameters setup 
取 值 单位 
ESC 工 况 〈 除 怠速 点 ) -- 
R123 -- 
3 MPa 
10 K 
0.4 MPa 
水 冷 〈 风 扇 辅 助 ) -- 
四 
eo 
OO 
| 
O 
B75 
J B75 加 
A/B/C25 A/B/C50 A75 A100 B100 C100 


50 100 | 150 200 250 300 | 350 
发 动机 的 功率 [kW] 

图 7 效率 模型 与 仿真 模型 的 对 标 结果 

Fig.7 Calibration of the efficiency model 
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虽然 无 法 获得 足够 的 试验 数据 验证 所 建立 的 


效率 模 


型 ， 但 是 如 果 所 建立 的 效率 模型 在 趋势 上 


完全 能 够 与 详细 的 仿真 模型 吻合 ， 有 潜力 代 蔡 详 


细 的 仿 
化 与 控 


真 模型 ， 用 于 ORC 余热 回收 系统 的 在 线 优 
制 。 


3.2 模型 应 用 验证 


将 
优化 控 
压力 ， 
合 特性 


所 建立 的 效率 应 用 到 图 4 所 示 ORC 系统 的 
制 架构 中 ， 在 线 优 化 决策 蒸发 压力 和 冷凝 
改善 ORC 系统 的 动态 特性 、 工 况 特性 和 耦 


o 


车 
工 况 ， 
需要 说 
是 有 限 
否则 起 
效 性 。 
容量 ， 


发 器 管 


辆 满载 ， 且 完全 采用 风扇 冷却 ， 在 HWFET 
对 所 建立 效率 模型 进行 应 用 效果 的 验证 。 
明 的 是 : HWFET 工 况 时 间 为 765s， 时 间 
的 ， 因 此 ， 需 要 控制 ORC 系统 的 热 容 量 ， 
动 时 间 太 长 ， 反 而 无 法 验证 控制 方法 的 有 
因此 ， 在 该 验证 过 程 中 ， 采 用 了 较 低 的 热 
为 50kg 不 锈 钢 当量 ; 在 该 条 件 下 ， 因 为 蒸 
壁 的 著 热 能 力 较 差 ， 对 排 气 扰动 的 缓冲 能 


力 不 足 ， 其 对 ORC 系统 的 蒸汽 状态 影响 较 大 ， 在 


这 样 比 


较 不 利 的 条 件 下 ， 验 证 所 提出 的 控制 架构 


和 效率 模型 ， 应 该 更 具有 普 适 性 。 验 证 结果 如 图 


8 所 示 


o 
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图 8 ORC 系统 在 HWFET 工 况 下 的 控制 效果 
Fig.8 Control effect of the ORC system during the HWFET driving cycle 


8， 对 ORC 余热 回收 系统 在 线 优 化 与 


控制 的 实现 过 程 阐 述 如 下 : 


1D) 


2) 


3) 


采用 低压 起 动 策略 ， 缩 短 ORC 系统 起 动 
时 间 。 

在 起 动 过 程 中 ， 冷 凝 压力 为 上 限 值 
0.8MPa， 在 台 架 运行 中 该 值 会 受到 工 质 
泵 入 口 压力 的 限制 。 该 决策 使 得 冷却 系 
统 不 工作 ， 一 方面 可 以 降低 冷却 水 泵 和 
冷却 风扇 的 功 耗 ， 另 一 方面 可 以 热量 的 
散失 ， 促 进 ORC 系统 快速 暖 机 和 起 动 。 
ORC 系统 在 起 动 、 冲 转 之 后 ， 进 入 到 做 
功 模 式 。 在 该 模式 下 ， 管 理 器 与 效率 模 
型 协同 工作 ， 完 成 三 大 功能 : 1 是 工 况 模 
式 的 管理 ， 尽 可 能 避免 ORC 系统 进入 到 
保护 模式 ; 2 是 动态 效率 的 优化 ， 在 一 段 
时 间 内 ， 控 制 ORC 系统 的 蒸发 压力 与 排 
气温 度 相 适 应 ， 提 高 ORC 系统 效率 ; 3 
僵 亏 看 合 关 系 的 平衡 ， 管 理 冷 凝 压 力 ， 
优化 调节 ORC 系统 发 电 与 冷却 系统 功 耗 


之 间 的 三 亏 关系 。 具 体 实现 结果 如 下 : 
a) ” 当 车 辆 突然 减速 ， 排 气 能 量 不 足 ， 
蒸汽 温度 下 降 ， 过 热度 也 会 下 降 ， 这 
时 ， 通 过 降低 蒸发 压力 ， 维 持 蒸汽 的 过 
热度 ， 避 免 其 在 短 时 间 内 跌 破 安全 限 值 
而 进入 到 保护 模式 。 当 然 ， 与 此 同时 ， 
控制 器 也 会 使 得 工 质 泵 的 流量 降低 ， 与 
蒸发 压力 协同 配合 ， 在 整个 循环 工 况 
中 ，ORC 系统 的 蒸汽 过 热度 一 直 控制 在 
10K 左右 的 水 平 。 
b) “ 当 车 辆 高 速 巡 航 ， 或 者 加 速 时 ， 排 
气 能 量 充足 ， 蔡 汽 过 热度 有 充分 的 保 
障 ， 在 这 种 情况 下 ， 蒸 发 压力 逐渐 升 
高 ， 保 持 ORC 系统 的 高 效 运行 。 

c) 针对 ORC 系统 发 电 与 冷却 系统 功 耗 
之 间 盔 亏 关系 的 平衡 问题 ， 当 车 辆 加 速 
时 ，ORC 系统 的 持续 吸 热量 和 散热 量 都 
较 大 ， 冷 却 系统 负荷 比较 大 ， 这 时 ， 通 
过 提高 冷凝 压力 的 控制 目标 ， 增 大 冷却 


构 
热度 
的 比 


系统 的 散热 温差 和 能 效 比 ， 降 低 冷 却 系 
统 功 耗 ， 当 车 辆 减速 或 者 巡航 时 ，ORC 
系统 的 持续 吸 热量 和 散热 量 较 小 ， 冷 却 
系统 的 功 耗 比较 小 ， 在 这 种 情况 下 ， 降 
低 冷 凝 压 力 的 控制 目标 ， 以 提高 ORC 系 
统 的 转化 效率 为 主 。 
总 体 而 言 ， 采 用 所 建立 的 效率 模型 和 控制 架 
ORC 系统 的 有 效 做 功 时 间 达 到 94%， 蒸 汽 过 
空 制 在 SK-15K 之 间 ， 散 热 功 耗 与 发 电功率 
值 维持 在 20% 以 下 ， 在 HWFET 工 况 中 的 节 


一 入 


效果 为 2.2% 一 2.8%。 


4 结论 


本 文 为 突破 ORC 余热 回收 系统 在 线 优化 控 和 


一 局 


的 瓶颈 ， 纤 合 考虑 ORC 系统 的 动态 特性 、 工 况 特 


性 和 
型 
架构 
构 的 


耦合 特性 ， 建 立 了 ORC 系统 的 效率 预测 模 
并 提出 了 基于 效率 预测 模型 的 在 线 优化 控制 
和 策略 ， 仿 真 结果 验证 了 效率 模型 和 控制 架 
可 效 性 ， 形 成 如 下 结论 : 
1) “柴油 机 ORC 余热 回收 系统 存在 强烈 的 动 
态 特 性 、 工 况 特性 和 耦合 特性 
于 稳 态 标定 的 优化 控制 律 在 动态 运行 过 
程 中 无 法 实施 ， 需 要 对 蒸发 压力 和 冷冻 
压力 进行 在 线 优 化 。 
2) ”以 ORC 系统 的 热力 学 第 一 定律 模型 为 基础 ， 抓 
住 蒸发 器 阶 跃 响应 的 宏观 特征 ， 提 出 了 有 效 换 
热量 的 动态 修正 模型 ， 以 及 工 况 模式 切换 的 预 
判 因 子 ， 并 采用 简化 机 理 和 数据 建 模 相 结合 
方法 ， 建 立 了 冷凝 器 散热 功 耗 的 经 验 模型 ， 旨 
成 建立 了 ORC 余热 回收 系统 的 效率 预测 模型 。 
3) ”所 建 效 率 预测 模型 在 仿真 平台 上 的 对 标 
误差 小 于 5%， 且 在 趋势 上 完全 能 够 与 详 
细 的 仿真 模型 吻合 ， 可 用 于 ORC 余热 回 
收 系统 的 在 线 优化 与 控制 。 
4) HWFET 道路 仿真 结果 从 应 用 效果 上 验证 
了 所 建 效率 预测 模型 的 有 效 性 ，ORC 系 
统 的 有 效 做 功 时 间 达 到 94%， 蒸 汽 过 热 
度 控制 在 5K-15K 之 间 ， 散 热 功 耗 与 发 
电功率 的 比值 维持 在 20% 以 下 。 
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